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ABSTRACT 

We present a comprehensive statistical analysis of Swift X-ray light-curves of Gamma- 
Ray Bursts (GRBs) collecting data from more than 650 GRBs discovered by Swift and 
other facilities. The unprecedented sample size allows us to constrain the rest-frame 
X-ray properties of GRBs from a statistical perspective, with particular reference to 
intrinsic time scales and the energetics of the different light-curve phases in a common 
rest-frame 0.3-30 keV energy band. Temporal variability episodes are also studied and 
their properties constrained. Two fundamental questions drive this effort: i) Does the 
X-ray emission retain any kind of “memory” of the prompt 7 -ray phase? ii) Where 
is the dividing line between long and short GRB X-ray properties? We show that 
short GRBs decay faster, are less luminous and less energetic than long GRBs in 
the X-rays, but are interestingly characterized by similar intrinsic absorption. We 
furthermore reveal the existence of a number of statistically significant relations that 
link the X-ray to prompt 7 -ray parameters in long GRBs; short GRBs are outliers 
of the majority of these 2-parameter relations. However and more importantly, we 
report on the existence of a universal 3-parameter scaling that links the X-ray and 
the 7 -ray energy to the prompt spectral peak energy of both long and short GRBs: 
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1 INTRODUCTION 

In 7 years of operation, Swift (Gehrels et al. 2004) has rev- 
olutionized our understanding of the X-ray emission that 
follows the prompt 7 -ray phase of Gamma-Ray Bursts 
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(GRBs): the X-ray afterglow. The standard afterglow the- 
ory (Meszaros & Rees 1997; Sari et al. 1998) predicts that 
X-ray emission arises from the interaction of a relativistic 
outflow with the ambient medium, leading to the forma- 
tion of a blast wave. In this context, short and long GRBs 
would naturally show similar afterglows (since the emis- 
sion would be sensitive to the energy budget of the rela- 
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tivistic outflow but otherwise keep no memory of the ori- 
gin of the outflow) if the properties of the environment of 
the two classes are also similar. This is difficult to recon- 
cile with the massive star (long GRBs, see e.g. Woosley 
1993) vs. compact binary (short GRBs see e.g. Paczynski 
1986) progenitor systems, which would instead suggest a 
wind density profile (i.e. oc r~ 2 ) around long GRBs and 
an ISM (i.e. oc r°) ambient density for short bursts. Con- 
trary to expectations, observations are often consistent with 
a constant density environment (ISM) also in the case of 
long GRBs (e.g. Racusin et al. 2009). Moreover, the stan- 
dard afterglow theory fails to explain the presence of long 
(~ 10 4 s) phases characterized by very mild decays (the 
so-called plateaus or shallow decay phases) in the X-ray 
light-curve of many GRBs (e.g. GRB 060729, Grupe et al. 
2007); it cannot account for abrupt drops of emission ob- 
served in some GRBs (e.g. GRB 070110, Troja et al. 2007, 
Lyons et al. 2010) and has serious difficulties explaining 
the X-ray flares (Chincarini et al. 2007, Falcone et al. 2007; 
Chincarini et al. 2010 and references therein). 

As a result, a number of alternative models have been 
proposed. They basically divide into two classes: accretion 
onto a newly born black hole which directly powers the ob- 
served X-ray light-curve (Kumar et al. 2008) or power from 
a rapidly rotating magnetar (see e.g. Metzger et al. 2011 and 
references therein). In sharp contrast to the standard after- 
glow theory, those models directly relate the properties of 
the observed X-ray light-curves to the GRB central engine. 

With this work we improve our understanding of the 
X-ray emission of long and short GRBs through a homo- 
geneous analysis of a sample of more than 650 GRBs ob- 
served by Swift-XRT (Burrows et al. 2005). We ask: what 
is the typical amount of energy released during the different 
X-ray light-curve phases? Does the X-ray emission retain 
any kind of “memory” of the prompt 7-ray phase? GRBs 
are traditionally classified into long and short according to 
their prompt 7-ray properties: do short GRBs show a dis- 
tinct behaviour in the X-rays as well? Is it possible to find 
a universal (i.e. common to long and short GRBs) scaling 
that involves prompt and X-ray properties? This set of still 
open questions constitutes the major reason to undertake 
the present investigation. 

Previous attempts mainly concentrated on observer 
frame properties and tried to understand to what extent 
the observations could be reconciled with the standard for- 
ward shock model (e.g. O’Brien et al. 2006; Butler 2007; 
Butler & Kocevski 2007; Willingale et al. 2007; Liang et al. 
2007, 2008; Evans et al. 2009; Racusin et al. 2009, 2011). 
This effort led to the identification of serious difficulties of 
the standard picture. We build on previous results and adopt 
here a different approach: instead of comparing observations 
with a particular physical model, we take advantage of the 
large sample size and look for correlations between the X- 
ray and 7-ray properties of GRBs any physical model will 
have to explain. We complement previous studies with: 

i. Homogeneous analysis of GRBs in a common rest frame 
energy band (0.3-30 keV). 

ii. Statistics and properties of the temporal variability su- 

perimposed on the smooth X-ray decay. 

iii. Comparative study of long vs. short GRB X-ray af- 
terglows. 


iv. Study of the prompt 7-ray vs. X-ray connection: no- 
tably, we report on the existence of a universal scaling in- 
volving prompt and X-ray parameters. 

Hereafter, we will refer to the X-ray signal recorded by 
Swift-X RT after a GRB trigger as “X-ray light-curve” (LC) 
and explicitly do not use the word “afterglow” to avoid con- 
fusion (“afterglow” refers to the standard interpretation). 
Uncertainties are given at 68% confidence level (c.l.) unless 
explicitly mentioned. Standard cosmological quantities have 
been adopted: Ho = 70kms^ 1 Mpc -1 , Da = 0.7, Dm = 0.3. 
The results from our analysis are publicly available. 1 


2 DATA ANALYSIS 

We select GRBs observed by Swift-XRT from the beginning 
of science operations in December 2004 through the end of 
2010. The starting sample includes 658 GRBs; 36 belong to 
the class of short GRBs. Following Margutti et al. (2011b), 
the short or long nature of each event is established using 
the combined information from the duration, hardness and 
spectral lag of its prompt 7-ray emission: a prompt 7-ray 
duration Too $ 2 s coupled to a hard 7-ray emission with 
photon index F < 1.5 and a negligible spectral lag are con- 
sidered indicative of a short GRB nature. 

For each GRB, the data reduction comprises four steps: 

• Extraction of count-rate light-curves (LC hereafter) in 
the 0.3-10 keV energy band. 

• Time-resolved spectral analysis. 

• Flux calibration in the observer frame 0.3-10 keV en- 
ergy band and luminosity calibration in the rest-frame 0.3-30 
keV energy range for the sub-sample of GRBs with known 
redshift. 

• LC fitting. 

The Swift-XRT data have been analyzed using the lat- 
est version of the HEASOFT package available at the time 
of the analysis (v. 6.10). For each GRB we started from 
calibrated event lists and sky images as distributed by the 
HEASARC archive 2 . The following analysis made extensive 
use of the XRTDAS software package. The additional au- 
tomated processing was performed via custom IDL scripts. 
Details on the procedure followed can be found in Margutti 
(2009). 3 Here we note that: 

(i) While agreeing on the major steps, the data extraction 
adopted here is slightly different from the methods presented 
in Evans et al. (2007) and Evans et al. (2009): however, af- 
ter comparing every single LC obtained with the two tech- 
niques, we find that the methods lead to consistent results. 

(ii) Our flux and luminosity LC calibration is based on 
a time-resolved spectral analysis which is able to capture 
the spectral evolution of the source with time. Uncertainties 
arising from the spectral analysis have also been propagated 
into the final flux and luminosity LCs (this is essential to 
compute the significance of positive temporal fluctuations 
superimposed on the smoothly decaying LC, see Sec. 2.1). 

1 A demo version of the website is currently available at 
http://www.grbtac.org/xrt_demo/GRB060312Afterglow.html 

2 http: / /heasarc. gsfc.nasa.gov/cgi-bin/W3Browse/swift. pi 

3 Retrievable from http://hdl.handlc.net/10281/7465 
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Figure 1. Cartoon showing the different LC models used in this 
work. Both axes use logarithmic units. 


(iii) The intrinsic neutral hydrogen absorbing column NHhg 
of each GRB was computed extracting a spectrum in the 
widest interval of time with no apparent spectral evolution. 
The best-fitting NHhg was used as frozen parameter in the 
time-resolved spectral analysis above. The Galactic contri- 
bution was frozen to the value in the direction of the burst 
as computed by Kalberla et al. (2005). 

Since we corrected for the Galactic and intrinsic ab- 
sorption, the final results are unabsorbed 0.3-10 keV (ob- 
server frame) flux LCs. For the sub-sample of GRBs with 
known redshift we furthermore extracted unabsorbed lumi- 
nosity LCs in the 0.3-30 keV (rest-frame) energy band (ex- 
trapolating the best-fitting power- law spectrum). We con- 
servatively use only s derived from optical spectroscopy and 
photometric redshifts for which potential sources of degen- 
eracy (e.g. dust extinction) can be ruled out with high con- 
fidence. The complete list of redshifts used is reported in 
Table B1 (175 GRBs in our sample have redshift). 


2.1 Light-curve fitting 

The X-ray LCs of GRBs for t > 60s consist of smoothly 
decaying power-laws or broken power-laws with X-ray 
flares superimposed. Here we concentrate on the underly- 
ing smooth component. 

We considered only GRBs whose statistics were good 



Figure 2. 0.3-10 keV unabsorbed flux LC of GRB 060312 with 
best-fitting model superimposed (blue solid line). The red crosses 
mark the data points identified as “excesses” during the semi- 
automatic fitting routine and subsequently removed from the fit 
to obtain the best-fitting power-law plus broken power-law model. 
Vertical dashed lines: best estimates of the break times obtained 
as described in the main text. 


Table 1. Number of GRBs per LC type. For each type, the num- 
ber of GRBs with complete LCs (C-GRBs, defined as promptly 
repointed GRBs f rep < 300 s whose fading was followed up down 
to a factor ~ 5 — 10 from the background limit) and with detected 
positive fluctuations (F-GRBs) with respect to the best- fitting 
smooth decay is also reported. GRBs are classified as either C or 
U-like, and either F or N-like. U-like bursts have truncated LCs, 
while N-like GRBs show no evidence for flares. This classification 
refers to the 0.3-10 keV (observer frame) LCs. 


Light- curve types 


0 

la 

lb 

Ila 

lib 

III 

Total Number of GRBs 

114 

89 

61 

133 

18 

22 

C-GRBs 

42 

61 

53 

121 

17 

22 

F-GRBs 

23 

16 

24 

48 

8 

10 


enough to allow us to extract a spectrum to convert their 
count-rate LCs into flux LCs (total of 437 GRBs out of 658). 
We first fitted the entire sample of flux LCs in the 0.3-10 keV 
(observer frame) energy band. We then focussed on the sub- 
sample of GRBs with redshift and performed a second fit 
using the LCs in the common 0.3-30 keV (rest-frame). 

Our semi-automatic fitting routine is based on the 
X 2 statistic and closely follows the procedure outlined in 
Margutti et al. (2011a). We fit the following models. Defin- 
ing: 

f (Ni ,ai,t) = Ni t~ ai (1) 

(where IVi and a\ are the normalization and the slope of 
the power-law, respectively) and 

g(N 2 ,a 2 ,ai 3 ,tb,s,t) = N 2 ^ — ^ + (c) ^ ( 2 ) 

(where N 2 is the normalization, a 2 and a 3 are the slopes 
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of the broken power-law; tb is the break time while s is the 
smoothing parameter), the fitting models can be written as: 


• Simple power-law (model 0): 

F = f(N uai ,t) (3) 

• Smoothed broken power-law (model la and lb for s < 0 
and s > 0, respectively): 

F = g(N 1 ,a 1 ,a 2 ,t b i,si,t) (4) 

• Smoothed broken power-law plus initial (model Ila) 
power-law decay: 

F = f(Ni,ai,t) + g(N 2 ,a 2 ,a 3 ,t b2 ,si,t) (5) 

or final (model lib) power-law decay: 

F = g(N 1 ,a 1 ,a 2 ,t b i,si,t) + f(N 2 ,a 3 ,t) (6) 

• Double smoothly joined broken power-laws (model III): 

F = g(N 1 ,a 1 ,a 2 ,tbi,si,t) + g(N 2 , a 3 , a 4 , t b 3 , s 2 ,t) (7) 


The model number follows from the number of break times. 
For model Ila (lib, III) the first (second) break time is 
defined as the time when the second (first) component out- 
shines the first (second) component. Figure 1 illustrates 
the different models, while Fig. 2 shows the result for 
GRB 060312 taken as an example: in this case the semi- 
automatic procedure identified two episodes of emission in 
excess of the smooth decay (model Ila). The number of 
GRBs for each LC type are listed in Table 1. We refer to 
the LCs as “type” 0, la, lb, Ila, lib and type III GRBs in 
the following. 

The best-fitting parameters together with their uncer- 
tainties and associated covariance matrix were then used to 
derive the 0.3-10 keV (observer-frame) fluence of the entire 
LC, from the Swift-XRT re-pointing time to the end of the 
observation. No temporal extrapolation was applied at this 
stage. Note that the contribution from significant positive 
fluctuations has not been included. The fluence of the dif- 
ferent LC phases as defined by the temporal breaks was also 
calculated (Fig. 1). Results are listed in Table B1 and Table 
B3. We then followed the very same procedure to fit the 0.3- 
30 keV (rest-frame) LCs: Table B5 reports the energetics in 
this energy range. 

The list of LC points flagged as “excesses” during the 
fitting procedure (e.g. red crosses of Fig. 2) constituted for 
each GRB the starting sample to look for significant posi- 
tive fluctuations with respect to the best fit. The information 
contained in the covariance matrix was used to derive the un- 
certainties associated with the residuals with respect to the 
best-fit (residuals were at this stage calculated on the entire 
LC). We first selected positive fluctuations with a minimum 
ltr significance. We furthermore require the positive fluctua- 
tions to show a rise plus decay structure-, this procedure au- 
tomatically excluded single data-points scattering from the 
best fit. GRBs showing (not showing) such structures were 
flagged as “F” (“N”) in Table B1-B5. GRB 060312 in Fig. 
2, with two rising and decaying structures superimposed on 
the smooth decay qualifies as “F”-event. The fluence (energy 
for known z) of those excesses was calculated by simply inte- 
grating the flux of each LC bin over the bin duration (after 
subtracting the contribution from the underlying smoothly 
decaying emission) . Errors were propagated accordingly and 


Short GRBs 

050724 051221A 051227 060313 061006 
061201 070714B 070724A 070809 071227 
080123 080503 080919 090510 090515 
090607 100117A 100816A 101219A 


Table 2. List of 19 short GRBs with complete LCs. GRBs 
with detected temporally extended emission are in boldface 
(Norris et al. 2011). 

can be used to quantify how significant is the presence of 
emission in addition to the smooth power-law decay in each 
GRB. The fluence (energy) of positive fluctuations detected 
during the different LC phases (e.g. steep-decay, plateau, 
normal-decay, etc.) was also derived and listed in Table B3 
(0.3-10 keV, observer frame) and Table B5 (0.3-30 keV, rest 
frame). For simplicity, in the following we will use the word 
“flare” to refer to statistically significant positive fluctua- 
tions detected on top of the smoothly decaying component, 
being however aware that different kinds of variability pos- 
sibly contribute to the detected “flaring activity” . 


3 LONG VS. SHORT GRBS PROPERTIES 

The analysis above reduces the X-ray LC of GRBs to a set 
of measured parameters: temporal slopes; break times; to- 
tal isotropic energy (fluence) and energy (fluence) associated 
with the different LC phases; flare energy (fluence); spectral 
photon index temporal evolution; intrinsic neutral hydrogen 
absorption. This constitutes an unprecedented set of infor- 
mation homogeneously obtained on the largest sample of 
GRBs to date and represents the natural sample to look for 
correlations among the parameters. 

While we report the best-fitting parameters of the en- 
tire sample (Tables B1-B5), in the following we restrict our 
analysis to GRBs with “complete” LCs, defined as those 
GRBs re-pointed by XRT at t rep < 300 s and for which we 
were able to follow the fading of the XRT flux down to a 
factor ~ 5 — 10 from the background limit (or, equivalently, 
t en d ^4x 10 5 s). These GRBs are flagged as “C” in Tables 
B1-B5. The number of “C”-like GRBs per LC morphologi- 
cal type is reported in Table 1. “U”-like GRBs have instead 
truncated LCs. Short GRBs with complete LCs are listed in 
Table 3. Short GRBs with extended emission are in boldface 
(see Norris et al. 2011). 

3.1 Median X-ray light-curve of long and short 
GRBs 

We select the sub-sample of C-like LCs of GRBs with 
redshift observed in the common rest-frame time interval 
10 2 — 10 5 s and 10 2 — 10 4 s for long and short GRBs, respec- 
tively. These criteria resulted in a sample of 79 long GRBs 
and 9 short GRBs (Fig. 3). We consider here 0.3-30 keV 
(rest-frame) luminosity LCs. 

We combined the best-fitting profiles of Sec. 2.1 to pro- 
duce a median luminosity LC of a long GRB. The result is 
shown in Fig. 3: the median luminosity LC roughly decays 
as ~ t _1 (with a milder decay ~ t~ 0 ' 9 for (0.5 < t < 4) ks 
and a steeper ~ t -1 ' 2 decay after ~ 4ks). With the possible 
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Figure 3. Best fitting profiles of 9 C-like short GRBs with rest- 
frame time coverage 10 2 — 10 4 s (colour lines) superimposed on the 
sample of long GRBs (grey lines). Black (blue) solid line: median 
Log(L) light-curve for short (long) GRBs; dashed lines mark the 
la dispersion. The LCs have been calibrated in the rest-frame 
0.3-30 keV energy band. 


Figure 4. 0.3-30 keV (rest-frame) X-ray energy as a function of 
redshift. Black (red) points: long (short) GRBs. Blue dashed line: 
empirically derived detectability threshold. For z > 2 we are not 
sensitive to GRBs with Ex iso 5) 10 51 erg. For z > 1 there is 
no evidence for an evolution of the upper bound of Bx,iao with 
redshift. 


exception of the shallower section, this is in rough agree- 
ment with the prediction of the standard afterglow theory 4 
(Meszaros & Rees 1997; Sari et al. 1998). 

The decay of the median LC is steeper for short (oc 
t -1 ' 3 ) than for long GRBs (oc f _1 ) in the rest-frame time 
interval 10 2 — 10 4 s; short GRB X-ray LCs are on average 
less luminous by a factor ~ 10 — 30 than long GRBs X-ray 
LCs. This conclusion holds also considering long GRBs in 
the same redshift bin. However, Fig. 3 clearly shows that the 
two samples slightly overlap (see also Gehrels et al. 2008). 
The steeper decay that characterizes short GRBs causes a 
progressive shift of their luminosity distribution towards the 
low end of the long GRB distribution. 

3.2 Energetics of long and short GRBs 

Table 3 reports the analysis of the parameter distributions 
derived from the LC fitting of Sect. 2. A complete list of 
symbols can be found in Appendix A. The observed Ex, iso 
distribution peaks at ~ 7 x 10 51 erg, typically represent- 
ing ~ 7% of the 1 — 10 4 keV (rest-frame) E^, i SO , Figure 4 
shows that we are not sensitive to the population of bursts 
with Ex, iso < 10 51 erg for z > 2 (so that the low energy 
tail of the Ex, iso distribution is currently under-sampled). 
This is likely a non-detectability zone, consequence of the 
Ex, iso oc E°)f so of Sec. 4. For z > 1 there is no evidence for an 
evolution of the upper bound of Ex, iso with redshift, which 
may suggest that ~ 10 53 erg is a physical boundary to the 
Ex, iso distribution (the record holder is GRB 080721 with 
Ex, iso ~ 10 53 erg). In this respect we note that maximum 
budget Emax ~ 10 52 erg 5 is predicted by magnetar mod- 

4 A steepening to ~ 1.5 — 2 is predicted after the jet-break time 
if the outflow is collimated into a jet (Rhoads 1999; Sari 1999). 

5 It is not a given that GRB 080721 violates this limit: Ex, iso 

represents the isotropic equivalent X-ray energy, an overestimate 

to the true value if the emission -as we believe- is beamed. 


els (Usov 1992). The same pattern is followed by the flare 
energy E^f: for z > 2 we are not sensitive to E^f < 10 5 ° 
erg. 

Observations suggest that the GRB X-ray LCs consist 
of two distinct phases (see e.g. Willingale et al. 2007): a first 
steep decay phase tightly connected to the prompt 7 -ray 
emission (Tagliaferri et al. 2005; Goad et al. 2006); and a 
second phase characterized by a flattening of the LC (with 
limited evidence for spectral evolution, see e.g. Liang et al. 
2007) followed by a “normal decay” phase. Type Ha GRBs 
(Fig. 1) clearly show the presence of both components, with 
energy E\,x = E\ and f? 2 ,x = E 2 + E 3 , respectively; for type 
la LCs, the lack of spectral evolution and the typically mild 
slope ai resembling 0:2 of type Ha lead us to identify E 2 ,x = 
E 1 +E 2 ; type lb GRBs show strong spectral evolution during 
the first LC segment: this together with the transition to a 
milder decay at tb, 1 lead us to define Ei,x = E\ and E 2 ,x = 
E 2 ; in type lib LCs the spectral and temporal properties 
of the first segment (with slope cti) strongly suggest that 
Swift-X RT caught the end of the prompt emission in the 
X-rays: we therefore define E\,x = E\ + E 2 and f? 2 ,x = 
E 3 ; the same is true for type III GRBs: in this case we 
define .Ei,x = E\ + E 2 , E 2 ,x = E 3 + E 4 . The two phases 
release comparable energy (see Table 3), with E\,x and E 2 ,x 
peaking at ~ 1.1 x 10 51 erg and ~ 4x 10 51 erg, respectively. 

In each distribution, short GRBs populate the low en- 
ergy tail: Ex, iso ~ 10 5 ° erg, which is approximately 2 orders 
of magnitude less than a typical long GRB. Figure 4 also 
shows that short GRBs are less energetic than long GRBs 
in the same redshift bin. A systematic difference between 
the -still poorly constrained- jet opening angles of long and 
short GRBs, with short GRBs being less collimated than 
long GRBs, could in principle mitigate this energy gap. If 
we compare the energy released during the two phases sep- 
arately (i.e. early steep decline vs. plateau plus subsequent 
decay), we find an indication that short GRBs are more en- 
ergetically deficient during the second phase then in the first 
phase, i.e. E^ rt /E l °" 9 ~ 0.014 and Eg^/E 1 ™ 9 ~ 0.054. 
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Table 3. Characteristic quantities describing the parameter distributions (number of elements (#), mean (m). median (M), standard 
deviation ( SD ), skewness ( SK )) and best-fitting values from a Gaussian fit (mean (//), standard deviation (<r), normalization (iV)). 
Fluences (S) are given in 10“ 6 erg cm -2 , energies (E) in 10 50 erg, fluxes (F) in 1 0 6 erg s 1 cm -2 , luminosities (L) in I O' 1 * ergs 
times (t) in s, hydrogen column densities (NH) in 10 22 cm -2 . Note that logarithmic (linear) units have been used in the upper (lower) 
half of the Table. We refer the reader to Appendix A for the exact definition of the parameters listed below. X-ray energies, luminosities 
and intrinsic times have been computed in the rest-frame 0.3-30 keV energy band. All the other X-ray quantities refer to the 0.3-10 
(observer frame) band. 



# 

m 

M 

SD 

SK 

M 


<J 


N 

Log(S-y) 

386 

0.17 

0.18 

0.61 

-0.07 

0.17 ±0.05 

0.60 

± 

0.04 

105.1 ± 7.7 

Log(Tgo) 

334 

1.58 

1.73 

0.65 

-1.18 

1.67 ±0.07 

0.62 

± 

0.06 

93.9 ±9.4 

Log(E 45 - 150 ) 

151 

2.06 

2.37 

0.93 

-1.12 

2.26 ±0.14 

1.02 

! ± 

: 0.1 

91.6 ± 9. 

Log(E 1 ^ so ) 

78 

2.88 

3.01 

0.91 

-0.89 

3.01 ±0.14 

0.85 

± 

0.12 

38.9 ±5.7 

Log(E pk) 

78 

2.64 

2.71 

0.52 

-0.76 

2.56 ±0.1 

0.47: 

±1 

3.069 

23.6 ±4.7 

Z/OpCZ/pk^so) 

85 

2.43 

2.51 

1.00 

-1.90 

2.51 ±0.19 

0.91 

± 

0.12 

61.0 ±9.3 

Log(T* f) 

138 

1.18 

1.33 

0.59 

-0.65 

1.26 ±0.11 

0.67 

± 

0.09 

53.8 ±7.4 

Log( NHhg) 

161 

21.6 

21.8 

1.21 

-3.54 

21.9 ±0.1 

0.62 

± 

0.09 

42.5 ±7.3 

Log(Sx) 

316 

-0.38 

-0.42 

0.62 

0.23 

-0.46 ±0.05 

0.57 

± 

0.04 

68.6 ±5.3 

Log(S^) 

115 

-0.81 

-0.76 

0.82 

-0.38 

-0.76 ±0.12 

0.89 

± 

0.12 

51.3 ±6.0 

Log(S lt x.) 

211 

-0.90 

-0.92 

0.80 

0.12 

-0.89 ± 0.10 

0.88 

± 

0.06 

97.2 ±7.9 

Log(S 2 ,x) 

316 

-0.58 

-0.61 

0.63 

0.14 

-0.57 ±0.07 

0.70 

± 

0.06 

102.2 ± 8.5 

Log(Sf\) 

62 

-0.75 

-0.69 

0.73 

-0.38 

- 


- 


- 

Log{st\) 

71 

-1.08 

-0.98 

0.90 

0.03 

-1.08 ± 0.17 

0.95 

± 

0.16 

36.7 ±5.6 

Log(E x , iso) 

126 

1.67 

1.84 

0.81 

-0.67 

1.82 ±0.08 

0.88 

± 

0.08 

31.1 ±2.5 

Log(E^) 

59 

1.25 

1.40 

0.97 

-0.64 

- 


- 


- 

Log{E li x ) 

86 

1.00 

1.04 

0.92 

-0.28 

1.10 ±0.11 

0.94 

± 

0.08 

40.7 ±3.6 

Log(E 2 ,x) 

126 

1.45 

1.63 

0.92 

-0.94 

1.63 ±0.11 

0.82 

± 

0.10 

63.2 ±7.0 

Log{El\) 

35 

1.13 

1.38 

1.00 

-0.78 

- 


- 


- 

Log{El\) 

38 

1.04 

1.14 

0.98 

0.01 

- 


- 


- 

Log(t\) 

155 

2.66 

2.56 

0.48 

1.06 

- 


- 


- 

Log(t{) 

155 

3.94 

3.93 

0.73 

0.19 

3.93 ±0.14 

0.8: 

±1 

3.12 

59.2 ±8.9 

L°g{F\) 

155 

-4.23 

-4.26 

0.83 

-0.14 

-4.18 ± 0.11 

0.89 

± 

0.12 

61.1 ±6.6 

Log(F { ) 

155 

-5.01 

-4.94 

0.80 

-0.12 

-4.99 ± 0.11 

0.80 

± 

0.08 

63.4 ±6.9 

Log{tf F ) 

62 

2.13 

2.03 

0.62 

1.27 

- 


- 


- 

Log(tf F ) 

62 

3.58 

3.48 

0.74 

0.29 

3.53 ±0.11 

0.75 

± 

0.09 

23.0 ±2.9 

Log(Li) 

62 

0.54 

0.73 

1.25 

-1.07 

0.79 ±0.19 

1.26 

± 

0.16 

46.6 ±5.5 

Log(Lf) 

62 

-0.47 

-0.19 

1.19 

-1.19 

0.04 ±0.33 

1.41 

± 

0.29 

38.2 ±6.8 

OL st 

213 

3.96 

3.56 

2.34 

3.95 

3.22 ±0.51 

2.34 

± 

0.35 

238.0 ±35.0 

^sh 

155 

-0.16 

0.18 

1.23 

-4.06 

0.27 ±0.14 

0.52 

± 

0.12 

61.0 ± 13.0 

a n 

204 

1.59 

1.38 

1.04 

8.30 

1.34 ±0.13 

0.49 

± 

0.11 

93.0 ± 19.0 


This argues against a beaming related explanation, since the 
jet opening angles of long and short GRBs are expected to 
be more similar at late than at early times. Short GRB light- 
curves decay faster than long GRBs in the X-rays, typically 
resulting in shorter observations ( tend ~ 10 4 s vs. tend ~ 10 5 
rest-frame): however, using the average L oc f -1,3 scaling 
above, we find E^ rt (t < 10 4 s) - E^ rt (t < 10 5 s). Thus 
the relatively lower measured energy of the later LC phase 
in short GRBs compared to long GRBs is not due to the 
shorter observations. 


3.3 Intrinsic neutral hydrogen absorption in long 
and short GRBs 

The distribution of the intrinsic neutral hydrogen columns 
NHhg is portrayed in Fig. 5, panel (c). The distribution 
of measured NHhg is found to have an average value 6 


Solar abundances are used to determine the best- fitting NHhg ■ 


m = 10 21,6 cm~ 2 . When modeled with a lognormal distri- 
bution, the best fitting mean and standard deviation are: 
/x = 21.9 ± 0.1, a = 0.6 ± 0.1, in agreement with the esti- 
mates by Campana et al. (2010, 2012) obtained on smaller 
samples. However, differently from Campana et al. (2010): 
(i) our sample contains a larger number of GRBs for which 
no evidence for intrinsic absorption was found (upper limits 
in Fig. 5); (ii) we find evidence for a larger population of 
highly absorbed (NHhg > 10 22 cm -2 ) GRBs at low redshift 
(z < 2). 

A trend for increasing NHhg with redshift is apparent in 
Fig. 5, panel (a)\ however, our sensitivity to small amounts 
of intrinsic absorption decreases with increasing redshift due 
to the fixed XRT band-pass, which explains the higher per- 
centage of upper limits in the 4 to 6 redshift interval, and 
is at least partially responsible for the observed trend. The 
sample is furthermore redshift selected, which implies a bias 
against highly extinguished GRBs. The severity of this bias 
is possibly redshift dependent, with a dependence which is 
difficult to quantify. 

Short GRBs map the low end of the NHhg distribution, 
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Figure 5. Left panel : (a): Intrinsic neutral hydrogen absorption 
vs. redshift for long and short GRBs (grey and red dots, respec- 
tively). 90% upper limits are marked with arrows. Median NHhg 
values in different redshift bins are indicated with filled trian- 
gles: for each bin, the error bars span the lcr NHhg dispersion. 
(b): redshift distribution of the sample of long (yellow) and short 
(red) GRBs. (c): Intrinsic neutral hydrogen distribution for long 
(black line) and short (red line) GRBs. The dashed histogram 
includes upper limits. In panels (b) and (c) a dashed line indi- 
cates the median value for the entire distributions ((z) = 1.82, 
(NHhg) = 10 21 ’ 8 cm -2 ). 


with an average absorption NHIhj 14 = 10 21 ' 4 cm -2 (mean 
of the logarithm of NHf(x; rt ). Their properties are however 
consistent with the intrinsic absorption of long GRBs in the 
same redshift bin. A KS-test comparing the NHhg distribu- 
tion of long and short GRBs with 0 < z < 1 reveals that 
there is no evidence for long GRBs to show higher NHhg 
when compared to short GRBs in the same redshift bin (KS 
probability of 34% ). A possibility is that we missed the pop- 
ulation of long GRBs with even higher NHhg at low a (see 
above). This would imply that GRBs with low optical ex- 
tinction but high NHhg are typical of the high-redshift uni- 
verse, only (Watson & Jakobsson 2012). We conclude that 
using the available data, caution must be therefore used to 
interpret the long GRB NHhg distribution as a proof of their 
association to star formation (Campana et al. 2010) unless 
this association is meant to be extended to short GRBs as 
well. 


4 PARAMETER CORRELATIONS 

Here we proceed to look for 2-parameter correlations in- 
volving both X-ray and 7-ray properties. From a blind anal- 
ysis we found 199 statistically significant correlations (out 
of 946). We focus on the physically interesting, correlations. 
The significance of each correlation is estimated using the 
R-index r, the Spearman rank p and Kendall coefficient K 
(Table 4). Only correlations for which the chance probability 
associated with at least one of the test statistics is < 1CP 3 
have been listed. As a general note: 


• No significant correlation is found to involve the rest 
frame Tgo, the intrinsic NHhg or the LC temporal slopes; 

• We re-scaled the LC temporal breaks th by the Tgo, 
adding new parameters to our list: y = tb/Tgo- However, 
the use of re-scaled properties did not improve any of our 
correlations and are therefore not included in the following 
discussion. 

The correlation coefficients and the best-fitting power- 
law parameters of each correlation are listed in Table 4. 
Our best-fitting procedure accounts for the sample variance 
(D’Agostini 2005). 

4-0.1 The link between the X-ray and prompt y-ray energy 

Fig. 6, (left panel) shows that Ex, i H o is directly linked to the 
isotropic energy released in 7-rays during the prompt emis- 
sion E 7 , i so . A similar result was found by Willingale et al. 
(2007) on a smaller sample of GRBs. Here we show for 
the first time how short GRBs compare to long GRBs: no- 
tably, all short GRBs but GRB 050724 are outliers of the 
long GRB relation, with Ex, iso for short GRBs a factor 
~ 50 below that for long GRBs and having large disper- 
sion (the Ex, iso distributions are almost distinct for long 
and short GRBs, as shown in Fig. 4). A clear exception 
is GRB 050724 (Barthelmy et al. 2005; Grupe et al. 2006) 
which had a bright and long-lived X-ray afterglow with a 
powerful late time re-brightening (Bernardini et al. 2011; 
Campana et al. 2006; Malesani et al. 2007). This difference 
may be understood in terms of a different radiative efficiency 
(where rj-/ = E 7 /(E 7 + Ex), being Ek the outflow kinetic 
energy) during the prompt emission between short GRBs 
and XRFs (X-ray Flashes, i.e. GRBs with E 7 ,i SO < 10 52 erg 
in Fig. 6): in this picture, r/l l hort > rjf RF . The two popu- 
lations are clearly distinct in terms of spectral peak energy 
during the prompt phase, with E^ ort > Ef k RFs (Fig. 6). 
This may suggest that rj-y anti-correlates with E p k: this is 
further investigated in Sec. 4.1.2. 

Short and long GRBs occupy different areas of the 
Ex, iso vs. Epk plane (Fig. 6, upper right panel) as well, 
demonstrating how the information from the X-ray LCs can 
be used to infer the GRB nature. Again, short GRBs fall 
below the long GRBs. 

4-0.2 The X-ray plateau and the prompt y-ray phase in 
long and short GRBs 

In the literature, the shallow decay (or “plateau”) is as- 
sociated with an LC phase generally characterized by a 
mild slope (and absence of spectral evolution in the X-rays; 
Liang et al. 2007): this can be identified in type Ha and type 
III light-curves. In type Ha (III) GRBs, this phase starts 
at ti = fbi (ti = thz) and ends at tt = tb2 (tf = tb3), 
with temporal slope 02 (0:3) and energy E2 (E3). Short 
GRBs are under-represented in the class of GRBs showing 
clear evidence of plateaus in the X-rays. Only 2 short GRBs 
(out of 19 with C-like LCs 7 , 10%) possibly have plateaus: 
GRB 051221A (T 90 = 1.4 s) and GRB070714B (T 90 = 3 

7 36 was the number of short GRBs in our starting sample; only 
19 of these have C-like LCs. 
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Table 4. From left to right: X and Y parameters to be correlated (the best fitting law reads: Log(Y) = q + mLog(X)); best-fitting 
parameters as obtained accounting for the sample variance (D’Agostini 2005): slope (m), normalization (q), intrinsic scatter (cr); errors 
are given at 95% c.l. The last six columns list the value of the Spearman rank p, Kendall coefficient K and R-index r statistics and 
relative chance probability p. For each parameter couple, values reported in the first line refer to the entire sample, while in the second 
line we restrict our analysis to the long GRB class. X-ray fluences, fluxes and observer frame times are computed in the 0.3-10 keV 
(observer frame) energy band; luminosities, energies and rest-frame times are computed in the 0.3-30 keV (rest-frame) energy band. 


X 

Y 

m 

<? 


a 


p 

p(p ) 

K 

P(K) 

r 

p(r) 

Ex,iso 

i?x L 

1.07 ± 0.03 

-4.1 ±72.2 

0.57 

± 

0.01 

0.70 

3 x 10- 10 

0.52 

2 x 10 -9 

0.79 

< 10 _1 ° 



1.10 ±0.07 

-5.9 ± 198 

0.58 

± 

0.01 

0.61 

3 x 10“ 7 

0.45 

6 x 10“ 7 

0.70 

lO" 9 

Ex., iso 

Li 

1.21 ± 0.06 

-15.6 ± 169 

0.85 

± 

0.01 

0.58 

6 x 10“ 7 

0.45 

2 x 10“ 7 

0.70 

3 x 10“ 10 



1.26 ±0.08 

-18.1 ± 206 

0.85 

± 

0.01 

0.55 

3 x 10“ 6 

0.43 

8 x 10“ 7 

0.69 

10" 9 

Ex, iso 

Li 

1.39 ±0.06 

-23.6 ± 172 

0.83 

± 

0.01 

0.63 

3 x 10“ 8 

0.49 

2 x 10^ 8 

0.75 

< 10 _1 ° 



1.37 ± 0.08 

-22.8 ±212 

0.84 

± 

0.02 

0.60 

3 x 10" 7 

0.46 

2 x lO" 7 

0.72 

lO-ro 

/RF 

Cf 

L t 

-1.23 ±0.03 

51.9 ±0.46 

0.77 

± 

0.01 

-0.80 

< lO” 10 

-0.60 

< 10 _1 ° 

-0.77 

< 10 _1 ° 



-1.24 ±0.03 

52.0 ±0.45 

0.73 

± 

0.01 

-0.82 

< lO" 10 

-0.62 

< lO" 10 

-0.78 

< 10 _1 ° 

Lf 

L/2,X 

0.52 ± 0.01 

26.8 ± 11.8 

0.47 

± 

0.00 

0.67 

2 x 10“ 9 

0.51 

5 x 10 -9 

0.80 

< 10 _1 ° 



0.50 ±0.00 

27.6 ± 10.8 

0.43 

± 

0.00 

0.65 

10" 8 

0.50 

2 x lO" 8 

0.81 

< lO" 10 

tq 

to 

X 

Lu.x 

0.42 ± 0.02 

29.1 ±43.2 

0.81 

± 

0.01 

0.42 

4 x 10“ 5 

0.29 

4 x 10“ 5 

0.45 

6 x 10“ 6 



- 

- 


- 


0.28 

6 x 10“ 3 

0.19 

0.06 

0.29 

6 x 10“ 3 

tf 

Ft 

-0.79 ±0.01 

-7.80 ±0.09 

0.45 

± 

0.00 

-0.69 

< lO" 10 

-0.50 

< 10" 10 

-0.74 

< 10 _1 ° 



-0.79 ±0.01 

-7.78 ±0.09 

0.45 

± 

0.00 

-0.69 

< 10" 10 

-0.50 

< 10~ 10 

-0.74 

< lO" 10 

E'y, iso 

Ex,i so 

0.79 ±0.01 

10.0 ±20.6 

0.39 

± 

0.00 

0.86 

< 10 _1 ° 

0.69 

< 10" 10 

0.88 

< 10 _1 ° 



0.67 ± 0.01 

16.5 ± 18.8 

0.29 

± 

0.00 

0.82 

< 10- 10 

0.63 

< lO" 10 

0.88 

< 10 _1 ° 

E ■ 

J — y ‘~/,lSO 


0.89 ± 0.05 

3.85 ± 148 

0.65 

± 

0.01 

0.64 

8 x 10~ 5 

0.48 

10" 4 

0.74 

3 x 10“ 6 



0.93 ±0.10 

1.83 ± 287 

0.62 

± 

0.02 

0.56 

10“ 3 

0.41 

10“ 3 

0.67 

6 x 10“ 5 

Ey t iso 

Lu,x 

0.67 ± 0.03 

15.9 ±91.3 

0.81 

± 

0.02 

0.71 

2 x 10“ 7 

0.56 

2 x 10“ 7 

0.64 

5 x 10“ 6 



0.56 ±0.04 

21.6 ± 126 

0.77 

± 

0.02 

0.62 

5 x 10“ 5 

0.48 

4 x 10“ 5 

0.54 

5 x 10“ 4 

E ■ 

J —''y,iso 

L 2 ,x 

0.92 ± 0.01 

2.96 ±33.5 

0.51 

± 

0.01 

0.76 

< 10“ 10 

0.59 

< 10~ 10 

0.85 

< 10 _1 ° 



0.74 ±0.01 

12.6 ±35.1 

0.44 

± 

0.00 

0.67 

10" 8 

0.51 

2 x lO" 8 

0.81 

< lO" 10 

E-y^iso 

L t 

1.06 ±0.08 

-8.86 ±227 

1.03 

± 

0.04 

0.54 

9 x 10“ 4 

0.41 

7 x 10“ 4 

0.70 

8 x 10“ 6 



1.05 ± 0.09 

-8.43 ± 260 

1.06 

± 

0.04 

0.50 

3 x 10“ 3 

0.37 

3 x 10“ 3 

0.68 

2 x 10“ 5 

L p k 

Ex., iso 

- 

- 


- 


- 

- 

- 

- 

- 

- 



0.98 ± 0.02 

49.5 ±0.15 

0.37 

± 

0.00 

0.63 

lO" 7 

0.46 

5 x 10- 7 

0.76 

< lO" 10 

L p k 

Ex, iso 

- 

- 


- 


- 

- 

- 






0.48 ±0.01 

27.0 ± 16.4 

0.44 

± 

0.00 

0.58 

2 x 10“ 7 

0.42 

3 x 10“ 7 

0.74 

< 10 _1 ° 

L p k 

L t 

- 

- 


- 


- 

- 

- 






0.86 ± 0.03 

2.35 ±87.6 

0.87 

± 

0.02 

0.50 

10“ 3 

0.39 

5 x 10 -4 

0.76 

7 x 10" 8 

L p k 

®2,X 

- 

- 


- 


- 

- 

- 






0.60 ±0.01 

20.3 ± 15.9 

0.43 

± 

0.00 

0.58 

lO" 7 

0.42 

3 x 10“ 7 

0.82 

< 10" 10 

S 7 

Sx 

0.77 ± 0.01 

-7.80 ±0.09 

0.45 

± 

0.00 

0.79 

< 10" 10 

0.59 

< lO” 10 

0.77 

< 10 _1 ° 



0.82 ± 0.00 

-1.58 ±0.10 

0.37 

± 

0.00 

0.78 

< 10 _1 ° 

0.58 

< 10~ 10 

0.78 

< 10" 10 

E ■ 

-*^7,130 

T llh 

L X 

0.71 ± 0.01 

8.53 ±30.9 

0.55 

± 

0.01 

0.66 

3 x 10 _1 ° 

0.49 

1 x 10 -9 

0.77 

< 10“ 10 


0.54 ±0.01 

17.5 ±29.6 

0.45 

± 

0.00 

0.55 

2 x 10“ 6 

0.40 

2 x 10“ 6 

0.70 

< 10 _1 ° 

E'y, iso 

T 10min 

L X 

0.93 ±0.01 

— 1.17 ± 0.01 

0.45 

± 

0.01 

0.87 

< 10 _1 ° 

0.67 

< 10 _1 ° 

0.88 

< 10 _1 ° 


0.78 ±0.01 

6.73 ±32.6 

0.40 

± 

0.00 

0.82 

< lO” 10 

0.63 

< lO" 10 

0.84 

< 10 _1 ° 


s, extended emission not included). The corresponding per- 
centage for long GRBs is instead ~ 37%. 


The luminosity at the end of the plateau phase Lf is 
directly related to the total energy released in the second 
LC phase i?2,x (Table 4): _E 2 ,x oc Lf ' 52 . It is interesting to 
note that of the two short GRBs, 070714B is a clear outlier, 
while 051221A is only barely consistent with the correlation. 
The peculiar GRB 060218 also shows a lower than expected 
L 2 ,x. Dainotti et al. (2008) first reported a correlation be- 
tween Lf and If for long GRBs. Here we confirm the cor- 
relation (with best-fitting Lf oc (t RF ) -1,2 Fig. 7) and show 
that the two short GRBs with clear evidence of plateau are 
not consistent with the same scaling. 


4-0.3 The link between the X-ray luminosity and the 
prompt 7 - ray energy release 

The X-ray luminosity of the LC, Lx(t RF ), correlates with 
the 7 -ray energy released during the prompt emission for 
any t RF between 100 s and 10 5 s. Here we arbitrarily se- 
lect two rest-frame times (Fig. 8 ) as an example. We find 
that the scatter of the correlation evolves with time, with 
the Lx(t RF ) vs. -E 7 ,iso being tighter at early times (see Fig. 
8 ). For this plot we require the GRBs to have been observed 
at those rest-frame times but relax the LC completeness 
requirement. No extrapolation of the observed LC is per- 
formed. 

At early times the LC luminosity tracks i? 7 i i so with 
limited dispersion around the best-fitting model Lx° mm oc 
E°'f so . Short GRBs tend to lie below the best-fitting law of 
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Figure 6. Correlations involving i? 7pso (1 — 10 4 keV, rest-frame), Ex, iso (0.3-30 keV, rest-frame) and the rest-frame prompt peak energy 
-Epk- Black dots (red triangles): long (short) GRBs. Dashed line: best- fitting power-law model for the entire short plus long GRBs sample. 
Dotted line: best-fitting model for the long GRB class, only. The coloured and hatched areas mark the 68% confidence region around 
the best fit. Short GRBs and outliers are named. 
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Figure 8. 0.3-30 keV (rest-frame) X-ray LC luminosity measured at 10 minutes (left panel) and llhr (right panel) rest-frame vs. £? 7pso 
Colour coding as in Fig. 6. The scatter of the correlation increases with time. 
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Figure 7. Luminosity at the end of the plateau phase, 0.3-30 keV 
(rest-frame) vs. end-time of the plateau. Colour coding as in Fig. 
6. Blue points are for type la LCs. 

the long GRB class. When compared to the same relation 
at much later times (11 hours) we find that: (i) the rela- 
tion is now more scattered, suggesting that the X-ray LCs 


are more directly linked to the prompt y-ray phase at early 
than at late times; (ii) while the relation is highly scattered, 
we note that all short GRBs of our sample lie below the 
long GRB relation: this is consistent with the steeper decay 
of the average short GRBs LC when compared to the long 
GRBs LC found in Sec 3.1. Our analysis therefore does not 
confirm the previous results from Nysewander et al. (2009), 
who found that short and long GRBs are consistent with the 
same Lx vs. i so scaling (note however that their _E 7pso is 
computed in a much narrower energy band). 

4-0.4 Observational biasses: temporal extrapolation 

The Swift re-pointing time t rep and end time of the observa- 
tions tend vary from GRB to GRB. Since Ex , iso is obtained 
by integrating the luminosity of each LC between t rep and 
tend, one may wonder what is the effect of using different 
integration times for different GRBs. This is quantified as 
follows. To estimate the amount of energy lost at the end 
of the observations, we extrapolated the best-fitting profile 
of each GRB up to 10 7 s (rest frame) and integrated the LC 
luminosity up until that time. Since GRBs may experience 
a jet break at late times (Racusin et al. 2009), this compu- 
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Figure 9. Impact of the temporal extrapolation of the observed 
LCs on the Ex iso vs. E 7 , aso correlation. Colour coding as in Fig. 
6. Light blue dots: Lx .iso values have been computed integrating 
the luminosity over a common rest-frame interval of time. 


tation may lead to an overestimate of the real energy lost. 
The amount of energy possibly lost at the beginning of the 
observations 8 is estimated by conservatively extrapolating 
backwards in time the best-fitting profile to the minimum 
rest-frame Swift re-pointing time of our sample, which is 
12.5 s. For GRBs with Tgo F > 12.5 s, we adopt Tgo F as the 
starting time for the integration to avoid extrapolating the 
luminosity to unrealistic values. This approach leads to an 
overestimate of the amount of energy lost before t rep for 
the large majority of GRBs, as can be seen comparing the 
extrapolated temporal profile we are adopting here to the 
Swift - BAT emission at the same rest-frame time (see e.g. the 
Swift Burst Analyser BAT plus XRT LCs of GRB 050724, 
Evans et al. 2010). The corrected Ex , iso is shown in Fig. 9 
(light blue points). Larger corrections (up to a factor ~ 9 for 
GRB 090510) are found to be applied to short GRBs. In spite 
of the very conservative approach we find that short GRBs 
are still either barely compatible or not consistent with the 
long GRB relation (as before), while the long GRBs relation 
is almost unaffected by this correction. 

We therefore conclude that in a logarithmic space the 
different rest-frame integration time used does not create 
or destroy correlations. The Ex , iso vs. E 7 , i HO correlation has 
been used here as an example: this result applies to all the 
relations presented in this paper. 


4.1 Multiparameter correlations 

We look here for correlations involving more than 2 param- 
eters (either from the X-rays or from the 7-rays). We first 
discuss the results from a Principal Component Analysis 
(PCA, Sect. 4.1.1) and then show the existence of a tight 3- 
parameter correlation directly linking Ex, iso, E 7 , i HO and E p k 
both in long and short GRBs (Sect. 4.1.2). 


8 Note that the sample of C-like GRBs we use to look for corre- 
lations was pre-selected requiring an observed time of re-pointing 
tobs < 300 s to minimize this effect. 


Table 5. The three most significant PCs (85% of the total vari- 
ance) projected upon E~ ti i so , E p \ r _, L p k, TJq 1 and Ex, iso- 



PCI 

40% 

PC2 

66% 

PC3 

85% 

K ■ 

J - /r y,zso 

-0.561 

0.141 

0.171 

Epk 

-0.448 

-0.300 

-0.630 

Epk 

-0.502 

-0.389 

- 

rpRF 
1 90 

-0.121 

0.794 

-0.512 

Ex, iso 

-0.466 

0.331 

0.588 


4.1.1 Principal Component Analysis 

The Principal Component Analysis (PCA) is a statistical 
technique designed to find patterns in data: it uses orthog- 
onal transformations to convert a set of possibly correlated 
variables into linearly uncorrelated (orthogonal) variables. 
Given a set of N events (GRBs in our case) described by M 
parameters, the PCA consists of the diagonalization of the 
covariance matrix: the eigenvenctors found are called prin- 
cipal components (PCs), while the eigenvalues consist of the 
variance associated with each PC (see Jolliffe 2002, for de- 
tails). We performed a standardized PCA as recommended 
when the parameters have widely different variances: each 
parameter P is replaced by P — ( Log[P ] — Log[P]) / a Log[p]- 
In this case the matrix to be diagonalized is not a covari- 
ance, but a correlation matrix. Calculations were performed 
using the statistical package R 9 . 

In Sect. 4 we showed that Ex , iso is the X-ray parameter 
that still keeps information from the prompt 7-ray energy 
release. We now investigate its relation to other prompt pa- 
rameters, specifically E p k, L p k, Tgo F (and E 7 , iso ), using the 
PCA. This set of parameters is measured simultaneously in 
44 GRBs. Table 5 reports the three most significant PCs 
(86% of the total variance) projected upon the original 5 
variables. Each variable roughly contributes with compara- 
ble weight to the first PC; the second PC is instead domi- 
nated by jfgQ F . The third PC relates Ex , iso with E p k. This 
result suggests that, while E p k, E p k, Ex, iso and E 7 , i SO are in 
some way physically related to one another (see Sect. 4.1.2), 
the duration of the 7-ray energy release represents an addi- 
tional degree of freedom to the system. 


4-1.2 A GRB universal scaling: Ex, iso, E 7 ,i SO and E p k 


We look for a 3-parameter correlation involving Ex, iso, E 7 ,i SO 
and -E p k- The three variables are found to be correlated (see 
Fig. 10) with the following best fitting law (obtained follow- 
ing the method by D ’Agostini 2005): 


Ex, iso 


1q(0.58±0.25) 


r p(1.00±0.06) 

7, iso 

p,(0.60±0.10) 

l ^pk J 


(8) 


where Ex, iso, E 7 , i so and E p k are in units of ergs and keV, 
respectively. The intrinsic scatter is <te x iao = 0.30 ± 0.03 
(la). We note that: 


• This relation expands on the well known E p k- E 7 ,i SO 


9 http://www.r-project.org/ 
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Figure 10. Three parameter correlation involving Ex, iso, E-f PSO and E p ^. Dashed line: best-fitting relation (Eq. 8): dot-dashed lines 
mark the 95% confidence area around the best-fitting law. Notably, long and short GRBs share the same scaling, with short and sub- 
energetic GRBs (like GRB 060218) occupying the same area of the plot. Inset: evolution of the efficiency parameter e = Ex iso/E~ i so as 
a function of the spectral peak energy of the prompt emission E p ^. A reference e oc E~//' e scaling has been marked with a blue dashed 
line. 


(Amati et al. 2008 and references therein) relation with the 
introduction of a third parameter (Ex, iso). 

• It combines information from the prompt and from the 
X-ray energy release which follows the prompt. While short 
GRBs are clear outliers of the E p k - E-, pso relation, they per- 
fectly fit into the Ex^so-E^^o-Ep^ relation: the importance 
of the 3-parameter relation is that it combines short and 
long GRBs on a common scaling. As a result, considering 
the entire short plus long GRB sample, the scatter is re- 
duced by the introduction of the third variable (the intrinsic 
scatter of the Amati relation of our sample of long and short 
GRBs is 0E pk = 0.37 to be compared to the intrinsic scatter 
of the 3-paramter relation on Epk which is <7B pk = 0.29). 
Restricting our analysis to long GRBs, we find <TE pk = 0.17 
both for the _E p k- E- u i Sa and for the 3-parameter correlation. 

• Short GRBs (like GRB 051221A) and sub-energetic 
GRBs (like GRB 060218) occupy the same region of the 
7?pk“ £7, iso- Ex , iso space. The same is true for the peculiar 
long GRB 060614, later re-classified as a possible short GRB 
(Gehrels et al. 2006). In general, GRBs seem to divide into 
two groups with “normal” long GRBs occupying the upper- 
right area; short and peculiar GRBs together with XRFs 
(e.g. 050416A, 060218, 081007, 060614 also have a spectral 
peak energy below 60 keV) share the same lower-left region 
of the plot. 

• The best-fitting slope of the Ex , iso vs. i BO relation 
of Fig. 6 reads: m = 0.79 db 0.01 (see Table 4). The sig- 
nificant departure of m from 1 implies the more energetic 
long GRBs to have a lower t = Ex, iso/ £7, iso oc l/E°f so , 
with short GRBs being outliers of this relation. Interest- 
ingly, equation 8 implies: e oc u 1 (iu , suggesting that the key 

Pk 

parameter determining the 7-ray to X-ray ratio is not E-, pso 
but the spectral peak energy E p k irrespective of the nature of 
the GRB (either long or short). This is clear from the inset 


of Fig. 10: the higher the prompt peak energy the lower the 
e (the GRB with _E p k ~ 10 4 keV is GRB 090510). We re- 
fer to Zhang et al. (2007) for a discussion of GRB radiative 
efficiencies derived from X-ray data. 

• The e(-Epk) scaling above can be interpreted as a phys- 
ical dependence of the radiative efficiency 777 on _E p k: r]i = 
E-//{E-/ + E k ) ~ E-f/Ex oc E-//E x as long as < Ek- 
This would imply % oc E p ff . See Fan et al. (2012) for a dis- 
cussion of this finding in the context of GRB photospheric 
models. Alternatively, a similar scaling could result for the 
long population if long GRBs with lower isotropic E lpao 
are less beamed than high energy GRBs during the prompt 
emission, but show otherwise similar beaming during the 
subsequent X-ray phase. In the first case, the e(7? p k) scal- 
ing would give direct information about the dissipative pro- 
cesses behind GRBs; in the second case, it would be an ob- 
servational effect, that nevertheless would provide valuable 
information about GRB jets and their opening angles. A 
complete and detailed discussion is beyond the scope of the 
present work and is provided by a companion paper (B12). 

• In Sec. 4.0.4 we showed that the different time intervals 
over which Ex , iso has been estimated do not severely affect 
the Ex, iso - F/7,iso“ Fpk ■ 


5 SUMMARY AND CONCLUSIONS 

We performed a comprehensive statistical analysis of Swift 
X-ray light-curves of 658 GRBs detected by XRT in the 
time period end of December 2004- end of December 2010. 
For the first time we present and analyse: (i) the properties 
of GRBs in a common rest-frame 0.3-30 keV energy band; 
(ii) we furthermore perform a comparative study of long 
and short GRBs; (iv) we cross-correlate the prompt 7-ray 
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properties and the X-ray LCs properties. We report below 
a summary of the major findings. 

From the spectral analysis of GRBs with redshift (Sec. 

3.3): 

1. We find evidence for high intrinsic neutral hydrogen 
absorption NHhg > 10 22 cm~ 2 even at z < 2. The average 
value for long GRBs is NHhg ~ 10 21 ' 9 cm~ 2 (mean value of 
the logarithm of the NHhg). 

2. Short GRBs map the low end of the distribution with 
mean NHhg ~ 10 21,4 cm -2 . However, there is no evidence 
for short GRBs to show a lower NHhg when compared to 
long GRBs in the same redshift bin. 

The analysis of 297 long GRBs with complete X-ray 
light-curves 10 reveals that: 

3. The average energy released in X-rays (0.3-30 keV, 
rest-frame) is -Ex, iso ~ 7 x 10 51 erg typically representing 
~ 7% of E 7 , i HO (Sec. 3.2). The two quantities are statisti- 
cally correlated: Ex, iso oc E°;f so . Also: Ex, iso oc E p k (Sec. 
4.0.1). 

4. The Ex, iso distribution does not extend beyond 10 53 
erg (Sec. 3.2) possibly suggesting the existence of a 
maximum available energy budget (the record holder is 
GRB 080721 with Ex, iso ~ 10 53 erg). Also, for z > 2 we 
are not sensitive to the population of GRBs with Ex, iso < 
10 51 erg, so that the low-energy tail of the Ex, iso distribution 
is currently under-sampled. 

5. The X-ray luminosity of the LCs at any rest frame 
time between 100 s and 10 5 s is found to correlate with 
E 7 ,i so (Sec 4.0.3): the scatter of this correlation increases 
with time which might suggest that early time X-rays are 
more tightly related to the prompt phase. 

In the case of short GRBs (19 have C-like LCs): 

6. The median luminosity light-curve of short GRBs (Sec. 
3.1) is a factor ~ 10 — 30 dimmer than long GRBs in the rest- 
frame time interval 10 2 — 10 4 s, has a steeper average decay 
(oc t -1,3 vs. oc i _1 ) and shows no evidence for clustering at 
late times (contrary to long GRBs). 

7. Short GRBs populate the low-energy tail of the Ex, i so 
distribution, with E^'iso ~ BtF^xTso an d an average 
Ex^iao* ~ 10 50 erg (Sec. 3.2). Short GRBs are more energy 
deficient during the second LC phase when compared to long 
GRBs. 

8. Short bursts are clear outliers of the Ex, iso — E 7 , i SO and 
Ex, iso — Epk relations established by the long population, 
with Ex) iso* a factor > 50 below expectations (Sec. 4.0.1). 
Short GRBs are also found to lie below the Lx lhr vs. E 7 , i SO 
relation established by the long class. 

9. Short GRBs are under-represented in the class of GRBs 
showing clear evidence of plateaus in the X-rays. Only 2 
GRBs out of 19 possibly have plateaus (10%). The corre- 
sponding percentage for long GRBs is instead ~ 37%. While 
the limited sample size does not allow us to draw firm con- 
clusions, we note that X-ray plateaus are more commonly 
detected in long GRBs (Sec. 4.0.2). 

10. The two short GRBs with X-ray LC plateaus in our 
sample are outliers of the Ef vs. If relation (Sec. 4.0.2). 

10 The total number of C-like GRBs is 316 (Table 1). 19 are short 
GRBs. 


Irrespective of the long or short GRB nature, we find no 
statistically significant correlation involving the rest frame 
prompt duration Tqq F , the intrinsic column density NHhg 
or the temporal slopes of the X-ray LCs (Sec. 4). The Tqq F 
basically accounts for the second strongest Principal Com- 
ponent (Sec. 4.1.1), suggesting that while E p k, E p k, Ex, iso 
and E 7 , i S o are related to one another, the 7-ray duration 
represents an additional degree of freedom to the system. 

We showed in Sec. 4.1.2 the existence of a 3-parameter 
correlation that links Ex, iso, E 7 , i SO and E p k: Ex ,iso OC 

/ rpl.OO / 77i0.60\ . 

(i) Short and long GRBs share the same scaling. 

(ii) This correlation implies ^' 13 ° oc E° jf which can be 
interpreted as r/ 7 oc E p k® (where r/ 7 is the radiative effi- 
ciency) . 

(iii) Standard long GRBs and short GRBs (together with 
peculiar GRBs and XRFs) occupy a different region of the 
Ex,iso - E 7 ,iso~E p k plane. 

The results from our analysis are publicly available. 11 
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APPENDIX A: GLOSSARY 

This section provides the list of symbols used. As a gen- 
eral note: X-ray energies (fluences) were computed from the 
time of the Swift-XKT repointing up until the end of the 
observation; no temporal extrapolation was performed. The 
values reported assume isotropic emission. X-ray fluences 
and fluxes are reported in the 0.3-10 keV (observer-frame) 
energy band; energies, luminosities and intrinsic time scales 
are computed in the 0.3-30 keV (rest-frame) band. 


• a n : temporal slope of the normal decay phase. Type la: 
Q n = ai\ type Ila: a n = < 23 ; type III: a n = 04 - 

• a s t: temporal slope of the steep decay phase. Type lb 
and Ila : a s t = an; type lib and III: a H t = 02 (see Fig. 1). 
The zero-time of the power-law decay is assumed to be the 
BAT trigger time (i. e. to = 0). 

• a™ 0 : temporal slope of the steep decay phase assuming 
to = Tgo . 

• a s h: temporal slope of the shallow decay (or plateau) 
phase. This corresponds to 02 and 0:3 for type Ila and type 
III light-curves, respectively. 

• Fj,: XRT 0.3-10 keV (observer frame) spectral photon 
index from this paper. 

• E-y, iso: isotropic equivalent energy released during the 
prompt emission in the rest-frame 1 — 10 4 keV energy band 
from Amati et al. (2008). 

• Epk: rest-frame peak energy of the vF v spectrum during 
the prompt 7 -ray emission from Amati et al. (2008). 

• Ff (Lf): flux (luminosity) at the end of the plateau (i.e. 
at t = tf). 

• Fi (Li): flux (luminosity) at the beginning of the plateau 
(i.e. at t = ti). 

• Lpk.iso: 1 — 10 4 keV (rest frame) isotropic peak lumi- 
nosity during the prompt emission from Nava et al. (2008). 

• JL 3 c h : luminosity at 11 hours rest-frame. 

• L 3 ? mln : luminosity at 10 min rest-frame. 

• NHtot: total neutral hydrogen column density. 

• NHhg: intrinsic neutral hydrogen column density at the 
redshift of the GRB. 

• 5i,x (Ei,x): fluence (energy) released during the first 
phase of the X-ray light-curve. Type lb and Ila: Ei,x = Ey, 
type lib and III: Ei,x — Ei + Eg. Fluences follow the same 
definition scheme. E\, Eg, Eg and E 4 has been defined fol- 
lowing Fig. 1 . 

• 52 , x (I? 2 ,x): fluence (energy) released during the second 
phase of the X-ray light-curve. Type la: E 2 ,x = E\ + Eg- 
type lb: Eg,x = Eg-, type Ila: Eg,x. = Eg + E 3 ; type lib: 
F/ 2 ,x = E 3 ; type III: l? 2 ,x = Eg + E 4 (see Fig. 1). Same 
definition scheme for fluences. 

• S-y (E 45 ~ 150 ): 15-150 keV (observer frame) fluence (en- 
ergy) released during the prompt emission as calculated by 
Sakamoto et al. (2011). 

• 5 X (E x , iso): X-ray fluence (energy). 

• 5x L (Ex L ): X-ray fluence (energy) associated to flares. 
For each GRB, the total fluence (energy) released in X-rays 
reads: 5£ L +5 X (E£ L +£ x ,i S o). 

• 5f ! x, 5 F ,x (Ef,x, Efjx): X-ray fluence (energy) of flares 
superimposed on the first and second light-curve phase. 

• tf, tf F , tj 90 : end time of the plateau phase: observer 
frame, rest frame and in T 90 units. This parameter corre- 
sponds to tb 2 and tb 3 for type Ila and type III light-curve, 
respectively. 

• ti, t; , t^ 90 : start time of the plateau phase: observer 
frame, rest frame and in T 90 units. This parameter corre- 
sponds to tb 1 and tb 2 for type Ila and type III light-curve, 
respectively. 

• At: plateau duration defined as tf — ti. 

• Tgo, 7c5q F : duration of the 15-150 keV prompt emission 
from Sakamoto et al. (2011), in the observer and in the rest- 
frame, respectively. 
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Table Bl. Best-fitting parameters of the 0.3-10 keV (observer frame) light-curves in de-absorbed flux units, as obtained following the 
procedure outlined in Sec. 2.1: GRB name, LC type (as defined in Sec. 2.1), redshift, power-law indices (a-i, a. 2 , < 23 , 0 : 4 ) and errors, 
break times (t^i, tb2> tb3 ) and errors, normalizations (Ni, N2) and errors, smoothing parameters (si, S2), prompt emission T 90 (we refer 
to Sakamoto et al. 2011 for GRBs detected before December 2009, and to the refined BAT GCNs otherwise), power-law index of the first 
segment when to = T 90 (c*^ 90 ) and error, x 2 , degrees of freedom, p-value. Normalizations are given in 10 — 10 erg cm -2 s — 1 , times in 
seconds. A redshift equal to 0 indicates that no reliable estimate of this parameter is available from the literature. For the other columns, 
—9 indicates that the value is absent (i.e. there is no such LC phase). Note on the LC type column: C-GRBs (i.e. GRBs with complete 
LCs) are defined as promptly repointed GRBs t rep < 300s whose fading was followed up to a factor 5-10 from the background limit; if 
this is not the case, the GRB is flagged as U-like (i.e. GRB with truncated LC). The flag F (N) indicates that flares have (have not) been 
detected. This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance. 


GRB 

Type 

z 

01 

&Oil 

OL2 

O' a.2 

cx. 3 

Oo 3 

Q4 

Oq4 

Log[tbi] 

Otbl 

Log[t b 2 ] 

Otb2 

Log[t b3 ] 

Otb3 

041223 

0UN 

0 

1.91 

0.44 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

050124 

0UN 

0 

1.44 

0.17 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

050126 

ICN 

1.29 

2.5 

0.51 

0.862 

0.25 

-9. 

-9. 

-9. 

-9. 

2.72 

0.4 

-9. 

-9. 

-9. 

-9. 

050128 

IUN 

0 

0.758 

0.12 

1.38 

0.14 

-9. 

-9. 

-9. 

-9. 

3.69 

0.38 

-9. 

-9. 

-9. 

-9. 

050219A 

ICN 

0 

3.68 

0.36 

0.779 

0.11 

-9. 

-9. 

-9. 

-9. 

2.4 

0.054 

-9. 

-9. 

-9. 

-9. 

050219B 

0UN 

0 

1.42 

0.044 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

050315 

IICN 

1.95 

-0.295 

0.18 

0.819 

0.037 

3.83 

0.18 

-9. 

-9. 

2.71 

0 

3.82 

0.14 

-9. 

-9. 


Table B2. Continued from Table Bl 


Ni 

ONI 

n 2 

ON2 

Si 

S2 

T90 

q T90 

O a T90 

x 2 

d.o.f. 

p-val 

1.45e7 

6.3e7 

-9. 

-9. 

-9. 

-9. 

109. 

1.91 

0.44 

8.39 

25. 

0.999 

4.72e4 

8.e4 

-9. 

-9. 

-9. 

-9. 

3.93 

1.44 

0.17 

31.6 

30. 

0.384 

0.0673 

0.11 

-9. 

-9. 

0.5 

-9. 

48. 

2.5 

0.51 

3.54 

10. 

0.966 

0.458 

0.49 

-9. 

-9. 

-0.5 

-9. 

28. 

0.758 

0.12 

69.8 

148. 

1 . 

0.457 

0.14 

-9. 

-9. 

1 . 

-9. 

23.8 

3.13 

0.31 

56.6 

79. 

0.973 

1.26e5 

4.8e4 

-9. 

-9. 

-9. 

-9. 

28.7 

1.42 

0.044 

81.9 

127. 

0.999 

0.179 

0.027 

2.15e9 

1.8e9 

-0.5 

-9. 

95.6 

-0.295 

0.18 

155. 

218. 

1 . 


Table B3. 0.3-10 keV (observer- frame) fluence table. From left to right: GRB name, LC type (as defined in Sec. 2.1), redshift, initial 
(' Tmin ) and final ( Tmax ) time of the observations, total fluence ( Sx iso ) with error, fluence of the different LC phases (Si, S 2 , S 3 , S 4 ) 
and errors, fluence of the flares in different parts of the LC ( S ^ L , S^ , S^ L ) and errors. Fluences are given in erg cm -2 . A 

redshift equal to 0 indicates that no reliable estimate of this parameter is available from the literature. A “—9” indicates that the LC 
does not show such phase and the value of that parameter is therefore absent. Finally, for the columns containing information from flares 
superimposed on the power-law decay, 0 indicates that no statistically significant positive fluctuation has been detected. This table is 
available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance. 


GRB 

Type 

z 

Tog[Tmin\ 

Log [Tmax] 

Sx, 

iso 

osx 


Si 

os 1 

s 2 

OS2 

S3 

OS3 

s 4 

OS4 

041223 

0UN 

0 

4.22 

4.45 

8.8e 

-8 

6.7e — 

9 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

050124 

OUN 

0 

4.05 

6.4 

1.6e 

- 7 

3.4e — 

8 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

050126 

ICN 

1.29 

2.12 

4.83 

3.7e 

-8 

6.8e - 

9 

1.5e — 8 

2.e — 9 

2.2e — 8 

6.3e - 9 

-9. 

-9. 

-9. 

-9. 

050128 

IUN 

0 

2.23 

4.84 

7.9e 

- 7 

2.4e — 

8 

4.6e - 7 

1.8e - 8 

3.3e — 7 

1.5e - 8 

-9. 

-9. 

-9. 

-9. 

05021971 

IUN 

0 

2.05 

6.2 

3.5e 

- 7 

1.4e - 

7 

4.1e — 8 

2.e — 9 

3.1e — 7 

1.4e - 7 

-9. 

-9. 

-9. 

-9. 

050219B 

0UN 

0 

3.5 

6.22 

9.4e 

- 7 

6.1e — 

8 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

050315 

IICN 

1.95 

1.92 

5.93 

l.le 

-6 

3.3e - 

8 

2.4e - 7 

1.4e - 8 

7.e — 8 

8.8e - 9 

8.e - 7 

2.8e - 8 

-9. 

-9. 


S[ L 

(J S1 FL 

qFL 

°2 

a S2 FL 

qFL 

^3 

a S3 FL 

qFL 

04 

a Si FL 

0 

0 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

0 

0 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

0 

0 

0 

0 

-9. 

-9. 

-9. 

-9. 

0 

0 

0 

0 

-9. 

-9. 

-9. 

-9. 

0 

0 

0 

0 

-9. 

-9. 

-9. 

-9. 

0 

0 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

0 

0 

0 

0 

0 

0 

-9. 

-9. 


Table B4. Continued from Table B3. 
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Table B5. 0.3-30 keV (rest-frame) energy table. From left to right: GRB name, LC type (as defined in Sec. 2.1), redshift, initial (T m j n ) 
and final (T marc ) time of the observations, total energy (Ex iso) with error, energy of the different LC phases (E\, E 2 , E 3 , E 4 ) and 
errors, energy of the flares in different parts of the LC (E ^ L , E^ 1, , E^ L , E^ L ) and errors. Energies are given in erg. A redshift equal 
to 0 indicates that no reliable estimate of this parameter is available from the literature. A “—9” indicates that the LC does not show 
such phase and the value of that parameter is therefore absent. Finally, for the columns containing information from flares superimposed 
on the power-law decay, 0 indicates that no statistically significant positive fluctuation has been detected. This table is available in its 
entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content. 


GRB 

Type 

Z 

Log[T m i n ] 

Lop [T maa; ] 

Ex, iso 

&EX 

E 1 

0\E1 

e 2 

&E2 

Es 

&E3 

Ei 

0\E4 

050126 

ICN 

1.29 

2.12 

4.83 

3.7e — 8 

6.8e - 9 

1.5e - 8 

2.e — 9 

2.2e — 8 

6.3e - 9 

- 9 . 

-9. 

-9. 

- 9 . 

050315 

IICN 

1.95 

1.92 

5.93 

l.le — 6 

3.3e - 8 

2.4e - 7 

1.4e - 8 

7.e — 8 

8.8e - 9 

8.e — 7 

2.8e — 8 

-9. 

- 9 . 

050318 

0UN 

1.44 

3.52 

5.65 

2.5e - 7 

l.le -8 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

-9. 

- 9 . 

050319 

IICN 

3.24 

1.98 

6.15 

5.e — 7 

5.e — 8 

2.5e - 8 

3.6e — 9 

1.7e — 7 

6.7e - 9 

3.e — 7 

4.9e - 8 

-9. 

- 9 . 

050401 

ICF 

2.9 

2.14 

5.9 

1.2e — 6 

6.2e — 8 

5.e — 7 

7.5e - 9 

7.3e - 7 

6.1e — 8 

-9. 

-9. 

-9. 

- 9 . 

050408 

0UN 

1.24 

3.41 

6.47 

5.4e - 7 

3.9e - 8 

-9. 

-9. 

-9. 

-9. 

-9. 

- 9 . 

- 9 . 

- 9 . 


Table B6. Continued from Table B5. 


e[ l 

&E1 fl 

E? L 

a E2 FL 

E el 

a E3 FL 

E^L 


0 

0 

0 

0 

-9. 

- 9 . 

-9. 

- 9 . 

0 

0 

0 

0 

0 

0 

-9. 

- 9 . 

0 

0 

-9. 

- 9 . 

-9. 

- 9 . 

-9. 

- 9 . 

0 

0 

0 

0 

0 

0 

-9. 

- 9 . 

0 

0 

0 

0 

-9. 

- 9 . 

-9. 

- 9 . 

0 

0 

-9. 

- 9 . 

-9. 

- 9 . 

-9. 

- 9 . 
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